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A nonlinear oscillator with an abruptly inhomogeneous restoring force driven by an uniform 
oscillating force exhibits stochastic properties under specific resonance conditions. This behaviour 
elucidates the elementary mechanism of the electron energization in the strong electromagnetic wave 
interaction with thin targets. 


1. INTRODUCTION 

Chaos in nonlinear dynamics plays a fundamental role 
in understanding the properties of various physical sys¬ 
tems pQ. A driven Duffing oscillator is an important 
example which demonstrates chaotic behaviour m- In 
particular, the problem of revealing the conditions for re¬ 
alization of chaotic regimes is of high importance for de¬ 
veloping the theoretical description of strong electromag¬ 
netic wave interaction with matter. This is connected to 
the mechanisms of charged particle energization in colli¬ 
sionless plasmas (this is also called ’anomalous heating’), 
in their interaction with high intensity electromagnetic 
waves. As an example we note here the well known elec¬ 
tron “vacuum heating” by laser radiation which occurs 
at the plasma-vacuum interface [3]. Such processes have 
been invoked as an elementary mechanism of collision¬ 
less absorption of electromagnetic waves and of electron 
heating during high power laser interaction with various 
targets (e.g. see j4]). Revealing of the elementary elec¬ 
tron heating mechanism is important for developing laser 
plasma sources of both fast electrons and ions. 

Among the ion acceleration mechanisms target normal 
sheath acceleration (TNSA) [5] has been the most inten¬ 
sively investigated. The TNSA mechanism crucially de¬ 
pends on the plasma electron heating. It implies that the 
accelerated ion energy is proportional to the mean elec¬ 
tron energy. The anomalous heating of electrons at the 
plasma-vacuum interface has been noted in Refs. Eld, 
where it was found that due to stochastic heating the en¬ 
ergy of the electrons recirculating around a small size 
target can be substantially higher than the energy of 
the electron quivering in the field of the electromagnetic 
wave. Taguchi et ah, El in analysing the results of the 
particle-in-cell simulations have called this regime “non¬ 
linear resonance absorption”. Finding the conditions and 
criteria of the realization of such a regime can have im¬ 
portant consequences, in particular, for optimizing the 
laser-target parameters in developing a laser ion accel¬ 
erator. Stochastic regimes in the electromagnetic wave 
interaction with a thin foil target have been found in 
theoretical papers M. We note here also the electron 


heating by a strong electromagnetic field due to the elec¬ 
tron rescattering off ions HM2. 

Below we analyze a model of a driven oscillator with a 
step-like nonlinearity in order to determine and elucidate 
the conditions when the stochastic vacuum heating takes 
place. The 3-dimensional model used bears the key fea¬ 
tures of the problem addressed in Refs. 0 inanu, but 
in a significantly simpler mathematical approximation, 
because a full description of the charged particle interac¬ 
tion with electromagnetic waves implies a consideration 
of a 7-dimensional dynamical system. In addition, such a 
mathematically simple system can lead to a broader un¬ 
derstanding of the underlying physical processes of elec¬ 
tromagnetic field interaction with charged particles, thus 
leading to developing novel theoretical tools. 

The paper is organized as it follows. In Section 2 
we formulate the mathematical model of the nonlinear 
driven oscillator and analyze the properties of its free and 
driven oscillations. We obtain the discrete mapping and 
find the threshold of stochastic regimes. The trajectory 
diffusion is described by the Kolmogorov-Fokker-Planck 
equation (KFPE) whose analytical solution is presented. 
Using the KFPE solution we derive the asymtotic depen¬ 
dence on time of the average electron momentum. Then, 
in Section 3, we present the result of the numerical inte¬ 
gration of the equations describing the trajectories of the 
driven oscillator with a step-like nonlinearity in nonrel- 
ativistic and relativistic limits, in regular and stochastic 
regimes. Here the time averaged electron kinetic energy 
is found numerically. At the end of the paper, in the 
Conclusion, we summarize the obtained results. 

2. DRIVEN OSCILLATOR WITH A STEP-LIKE 
NONLINEARITY 

We consider the system of ordinary differential equa¬ 
tions 

p + £ p sign(x) = a cost, (1) 

P 

x = -T IS' 

(i+p 2 ) 1/2 


( 2 ) 
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Here the sign function sign(x) = 1 for x > 0 and 
sign(x) = —1 for x < 0, and the dots denote differen¬ 
tiation with respect to time, t. The driver force on the 
r. h. s. of Eq. 0 represents the electromagnetic wave 
within the framework of the dipole approximation. 

In the case of the electron circulating around a thin 
foil target, p and x are the electron momentum and co¬ 
ordinate normalized by m e c and c/cj, respectively, where 
m e is the electron mass, uj is the driver field frequency, 
c is the light speed in vacuum. In Eqs. 00 the time 
variable is normalized by a; -1 . The parameter 


27rnoe 2 /o 

m e ujc 


( 3 ) 


introduced in Ref. m, is proportional to the charge 
separation electric field, E cs = 27reno/o, normalized by 
m e cjc/e, for a thin foil with the charge surface density 
equal to en^lo. The driver field normalized amplitude on 
the r. h. s. of Eq. 0 equals 


eE 0 

a = -. 

m e ujc 


( 4 ) 


If the driver field vanishes, a = 0, the system of equations 
00 describes free oscillations of a relativistic oscillator 
with an abruptly changing restoring force. The energy 
integral of Eqs. 00 is given by 

(l +p 2 ) ^ + £p\x\ = £. (5) 


Using this integral, we find the momentum and coordi¬ 
nate dependence in time for the case of free oscillations: 


P = Pm + sign(x)e p t, 


(6) 


* = xo + signer) g [l + (Pm (7) 


with the constant £ on the r. h. s. of Eq. 0 equal to 


The solution of Eqs. 01 can be written in terms 
of piecewise-smooth functions. Let the initial condition 
at t = to be x(to) = xo, and p(to) = pq. Let tk be 
consecutive zeros of the function x(£), k = 1, 2, 3,... (it 
is possible that t\ = to, but not required). 

The solution to p(t) in the interval t G [t n ,t n +i] is 


pit) =p 0 +£pS 0 (2 yy-i ) k t k + 1 0 + (-i ) n+i tj 

+a(sint — sinto)(ll) 


k=l 


where 


So={ 


' sign(a;o), + 0, 

sign(po), x o = 0,Po^0, 

sign(cos(t 0 )), x 0 =p 0 = 0, cos (t 0 ) 7 ^ 0, 
k -sign(sin(t 0 )), x 0 = po = cos (t 0 ) = 0. 


Then for x(t) we have 


t 

x(t) = j 

to 


Pit) 


(i +p{tyy/* 


dt. 


( 12 ) 


(13) 


The consecutive zeros of the function x(t) are deter¬ 
mined by the implicit relationships 

tk + l 

J UTj$W2 dt -°’ k = 0 ’ 1 ’ 2 ’---- < 14) 

tk 

Loops are possible, if there exist integers n > 2 and 
m > 1, such that the momentum at t n is the same as at 
ti, p n ~ Pi = 0, and the phase of the oscillating compo¬ 
nent of the force is the same, t n — t\ = 2irm. This leads 
to the condition 

n 

2^(-l) /c t /c -ti-h(-l) n+1 t n = 0, t n = t\ -\-2nm. (15) 
k =2 


£— [l + (pm ~ £p t] 2 ) ^ E s p \xq\. (8) 


It is determined by the initial values of the momentum 
p m and coordinate xq. 

Using relationships ([6| and 0 we find the period of 
nonlinear oscillations 



£p 


( 9 ) 


depending on their amplitude p m . 

Now we consider the case of driven oscillations, when 
a / 0, i.e. the electron moves in the superposition of the 
the static inhomogeneous electric field, — £ p sign(x), and 
oscillating homogeneous field, a cost. This is described 
by the system of equations (l]2) with the Hamiltonian 


l~L(x,p,t) = (1 + p 2 ) 1 / 2 + e p \x\ — ax cost. (10) 


While the electron trajectory crosses the x = 0 plane, 
the sign of the static electric field changes abruptly. This 
can be considered as a “collision”, during which the os¬ 
cillating driver electric field produces work changing the 
electron energy. Depending on the phase, when the “col¬ 
lision” happens, the work can be either positive or neg¬ 
ative, resulting in the electron energy increasing or de¬ 
creasing. 

A characteristic electron momentum change during the 
x = 0 plane crossing can be estimated to be of the order 
of 


A p — a. (16) 

The enhancement of the energy transfer from the 
driver electric field to the electron is expected to occur 
when a resonant like condition takes place, i. e. when the 
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periods of the motion of the electron in subsequent peri¬ 
ods of the free oscillations with momenta p and p + A p 
are co-measurable with the period of the driver electric 
field, 27r. This yields 


P 


J 


and 


P + Ap _ 7T 

£r> 2 


(17) 


(18) 


with j = 1, 2, 3,... and k = 2,3,4,.... From these expres¬ 
sions it follows that the resonance condition is 


7 r(k — j) 
a = -- —-£„ 


(19) 


This gives a relationship between the parameter e p and 
the normalized driver amplitude, a, which can be rewrit¬ 
ten in terms of the static and oscillating electric fields as 
E cs ~ Eq. Here and further it is assumed that k — j = 1. 


We note here that the condition (19) is equivalent to the 


criterion of thin slab relativistic transparency mm as 
well as determining the optimal parameters for ion accel¬ 
eration by the laser light radiation pressure mm- 
Using Eqs. (16-19) we can find the electron mo¬ 


mentum time dependence in the regimes of regular and 
stochastic acceleration. 

In the regular acceleration regime, the electron energy 
increases with A p = a during each half-period of free 
oscillations. Taking into account the free oscillation pe¬ 
riod dependence on the electron momentum © we obtain 
that the maximum of the momentum grows in proportion 
to the square root of time, 


p m = ( aSpt/2 ) 1 / 2 = a (t/ 7r) 1/2 . 


( 20 ) 


In order to describe the stochastic acceleration let us 
assume that the trajectory crosses the x = 0 plane at the 
instants of time equal to t n . The momentum change is 


Pn-\- 1/2 — Pn £p(tn+l/2 Ui) T CL sin(t n ), 


( 21 ) 


Pn+1 — Pn+1/2 + £p(tn+l ~ £ 71 + 1 / 2 ) + a s i n (Ui+l/2)- (22) 

In the limit of large momentum 

\Pn\/d£p^>l (23) 

the time interval between two subsequent crossings of the 
x = 0 plane is equal to 


^n+l tn 


4|Pn+l| 


(24) 


Here we have neglected the contribution of the order of 
a. Introducing the phase 0 n = t n we rewrite expressions 


(21-24) in the form of the Poncare map: 

Pn +1 =Pn + 2asin(0 n ), 


(25) 
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FIG. 1. Phase portrait (<p,p) of the mapping ( |25|26 ) for e v — 
1, po = 1 and 0o = 7 t/ 3. a) The stochastic layer formation 
for a — 0.25. b) The stochastic sea a = 0.3. 


<t>n+ 1 = 4>n + 4 ^ w+1 ^ (mod 2tt), (26) 


which is also known as the standard map, corresponding 
to the model describing stochastic Fermi acceleration [21 , 
If22] as a material point bouncing between two oscillating 
walls, e. g. see Hang. As is well known a Jacobian of 
the mapping is equal to unity, 


(dp n +l,d4>n+l) 

(dp n ,dcj) n ) 


= 1. 


(27) 


i.e. the mapping conserves the phase volume. Due to 


the condition of the driver amplitude smallness (23) the 


momentum change during subsequent crossing of the x = 
0 plane is relatively small while the phase <p can change 
significantly. Using this fact we can obtain from Eqs. 


(25) and (26) that 


0n+1 ~ + 


4|Pn| . 8|p n |a . 

-1-sin (<£„). 


(28) 


This expression yields for the phase volume stretching 
along the 0 coordinate 


\K\ 


1 


S<t>n 


- 1 


8a\p n 


COS (0 n ). 


(29) 


Since in the case of ’’highest probability” | cos(0 n )| ~ 1, 
the criterion of the stochastic regime onset is \K\ > 1 or 


SclIpyi 


> 1. 


(30) 


In Fi g. [T| wo show the phase portraits (0,p) of the map¬ 
ping ( 25|26 ) at the stochasticity threshold (a), and in the 
stochastic regime (b), when the particle energy can grow. 
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Eq. (25) yields for the averaged squared deviation of 
the momentum change during one cycle 


In the limit of small but finite momentum this expression 
yields 


(A p) 2 = (p n+ 1 - p n ) 2 = 4a 2 sin 2 (<^„) = 2a 2 . (31) 

According to Eq. 0 the cycle duration is equal to r = 
4\p\/e p . This gives 


(A p ) 2 = 


a 2 s p t 

2\p\ 


(32) 


which shows that the momentum dependence on time has 
a character of diffusion with the diffusion coefficient 


D 


ApAp 


pp 


a £ p 

2\p\ ' 


(33) 


Introducing the distribution function /(£,p), which 
obeys the KFPE, we write this equation as 

dtf = ( D pp d p f ), (34) 

or 

a/=%>„(y„/). (35) 

The solution of the KFPE for the initial condition 
f(t = 0,p) = no5(p), where S(p) is the Dirac delta func¬ 
tion, can be easily found. It reads 


f = n 0 


3 1 / 3 


where 


(8xi) 1/3 r(l/3) 


a Erj 


exp 


-S <»> 


(37) 


and T(x) is the Euler gamma function [25] . 

Using Eq. (36) we can find the time dependence of the 
average electron momentum, 


P= j Pf(p)dp■ 


It is given by 

3 2 / 3 p(2/3) vj /3 _ 3 2/3 r(2/3) 2 .\i /3 

p 2r(i/3) ^ ^ 2 4 / 3 r(i/3) ^ pt ^ ’ 

with the momentum dependence on time p ~ t 1 / 3 . 
The average kinetic energy, 

/*+oo 


m e c z 


/ too 

[(1 +p 2 ) 1/2 - 1 }f(p)dp 

-OO 


£, 


m e c z 


G\ 


1 + 


( 1 

1/6,1/2,5/6 \ 

\^324Ot) 2 

-1/3,0,0,1/3,1/2 ) 


48 x 3 1 / 6 7r 3 (>rf) 1 / 3 


(38) 


(39) 


(40) 


can be written in terms of the Meijer G - function 


( 41 ) 


4 „ r(i/3)r(i/6)r (5 /6) /3 = i ei , 1{xtf ^ 

m e c 2 47r 2 

(42) 

When t -MX) we have 


— ~ 2 2/3 7T 3/2 r(l/6) t \i/3 = 1 051(xf) 1/3 

3 4 / 3 r(i/3) 2 r(2/3) 2 ^ ’ • ( ) • 

_ (43) 

t 1 / 3 , has been 


m e c 


The energy dependence on time, £ e 
noted in Ref. [§]. 

The above derived standard map given by Eqs. 
( 25|26 ), the KFPE equation (35), and the KFPE solu¬ 
tion (36), showing the estimations for the stochasticity 


threshold (30) and the heating rate (39), have been ob¬ 


tained by using certain assumptions, which correspond 


to the conditions imposed by Eqs. (16, 23, 28). Since, as 


is well known, the properties and behaviour of dynami¬ 
cal systems in the stochastic regimes are very sensitive to 
the parameters, below we present the results of numeri¬ 
cal integration of the equations of motion in classical and 
relativistic limits. They provide additional information 
on the properties of a driven oscillator with a step-like 
nonlinearity. 


3. RESULTS OF NUMERICAL SOLUTIONS OF 
THE EQUATIONS DESCRIBING DRIVEN 
OSCILLATIONS 

In order to analyse in more detail the chaotic v.s. reg¬ 
ular dynamics of the driven nonlinear oscillator we have 
integrated Eqs. @§ numerically for different param¬ 
eters. In order to regularize the singularity on the left 
hand side of Eq. 0 we use instead of sign(x) the function 
Tanh(x//) with the width l substantially smaller than the 
typical particle displacement. In dimensional variables 

1 /o 

this requires l eE/m e cj 2 e , where uj pe = (47rne 2 /m e ) 
is the Langmuir frequency. 

The oscillator dynamics exhibits principally different 
behaviour depending on the relative values of the pa¬ 
rameter £ p and driver amplitude a. When the normalized 
driver amplitude is substantially less than unity, a « 1, 
we have the non-relativistic limit with nonlinear effects 
due to the abrupt dependence of the restoring force on 
the coordinate. In opposite limit, when a 1, in ad¬ 
dition there is a nonlinearity due to the relativistic rela¬ 
tionship between the velocity and momentum §. 


Non-relativistic regime 

In the classical limit the behaviour of the driven non¬ 
linear oscillations is determined by one dimensionless pa- 
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rameter, b = a/s p , which can be written as 


b = 


E 0 

2tt enl 


(44) 


This is equal to the ratio of the driver electric field to the 
electric field formed by a thin foil with the electric charge 
surface density 2iTenl. The system of equations 0 and 
Q takes the form 

x + sign (x) = bcos(t). (45) 


Here the x variable is normalized by c/ujs p = 
m e c 2 /2'Ke 2 nl = d 2 /1 with d e = c/uo pe being the collision¬ 
less skin depth. We shall also use the notation p = x. 

In the case of free oscillations when the r. h. s. of Eq. 
(45) vanishes, i. e. b = 0, the energy integral 0 takes 
the form 


d + \ x \=£ (46) 

with the constant £ = p%/2 + \xo\. The period of free 
oscillations is given by Eq. (J9|, which in the notations 
used in this subsection is T = 4 p m . 

The resonance condition ( p~9] ) in the non-relativistic 
limit yields 


b = 7r/2. 


(47) 


If the driver amplitude is substantially smaller than 
unity, for the initial conditions po = 0 and xo = 0 the 
trajectory is confined in the vicinity of the equilibrium 
x = 0 within the region \x\ <C l. 

For the driver amplitude approaching the threshold 
© the oscillations become nonlinear with the excur¬ 
sion exceeding the width l as seen in the case shown in 
Fig. m when the chosen parameters are close, but below 
the threshold (47). The dependence of the momentum 
and coordinate on time in Figs. [5] (a) and (c) correspond 
to nonlinear oscillations with slowly varying period and 
amplitude. In the phase plane, Fig. [2] (b), it is seen that 
the trajectory is confined within finite domain. In Fig. 
[2] (d) we plot the Poincare section showing the particle 
positions in the phase plane (x,p) at the discrete time 
with the time step equal to the period of the driving 
force, 27r. The loop trajectory is slightly broadened due 
to approaching the stochastic regime. 

There is a difference between the Poincare sections pre¬ 
sented here, and below in Figs. |(d),@(b,d,f),|5](d), 
i (d), 1 (a), [TT] (b), and in [l2] (b), and those Poincare 
sections which are shown in Figs. |T] (a) and (b). The 
driven oscillator dynamics is characterized by two time 
scales: the period of the driver force and the time in¬ 
terval between subsequent crossing by the trajectory the 
plane at x = 0. While in Fig. [l] and in Fig. [8] (b) each 
point in the plane (</>,p) corresponds to the particle mo¬ 
mentum and the phase calculated for the time when the 



FIG. 2. Driven oscillations below the stochastic regime 
threshold in the case of b — 0.775, l = 0.000125, po = 0, 
and xo = 0. a) The momentum p v.s. time, b) The phase 
plane (x,p) for t G [0,10 X 2n\. c) The configuration space 
(t,x,p) for t G [0,10 x 2tv\. d) The Poincare section showing 
the particle positions in the phase plane (x,p) at the discrete 
time with the time step equal to the period of the driving 
force, 2tv. 


trajectory crosses of the plane at x = 0, the points in the 
phase planes (x,p) presented in Figs. i(d),i(b,d,f),[5] 
(<*),© (d), [TT] (b), and in [l2] (b) show the particle coordi¬ 
nate and momentum at the subsequent instants of time 
separated by the driver force period. 

For the driver amplitude slightly above the stochastic- 
ity threshold being equal to b = 0.795 the oscillations 
become apparently stochastic as clearly seen in Fig. [3] 
During the initial interval of time f G [0, 50] the oscilla¬ 
tions have a relatively low, gradually growing amplitude 
as seen in Figs. [ 3 ] (a)-(c). Then it abruptly, during 
several periods, increases up to the level approximately 5 
times higher than the driver amplitude. The irregular os¬ 
cillations can roughly be represented as a superposition 
of two modes with higher and lower amplitude. There 
is also seen an intermittency when for some periods of 
time the oscillation amplitude becomes abruptly small 
and then increases again. The Poincare section in Fig. [3] 
(d) demonstrates the trajectory broadening and appear¬ 
ance of islands as is typical for stochastic regimes pH. 

Further increase of the driver amplitude results in de¬ 
veloped stochastic oscillations illustrated in Fig. [4] If the 
driver amplitude is equal to 1, the corresponding phase 
plane and Poincare section (not shown here) are similar 
to ones presented in Fig. [3] If the driver amplitude, 5, is 
equal to 5 the trajectory in the phase plane, Fig. [ 3 ] (a), is 
confined in the gradually broadening stripes. When the 
driver amplitude, 6, is equal to 9 the trajectory tightly 
fills the confinement domain in the phase plane (see Fig. 
[ 4 ] (c)) showing ergodicity. The Poincare sections in Figs. 
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FIG. 3. Driven oscillations for the driver amplitude higher 
than the stochasticity threshold for b — 0.795, l — 0.0125, 
po = 0, and xo = 0. a) The momentum p v.s. time, b) The 
coordinate x v.s. time, c) The phase plane (x,p). d) The 
Poincare section showing the particle positions in the phase 
plane (x,p) at the discrete time with the time step equal to 
the period of the driving force, 27r. 

[4] (b) and (d) demonstrate the formation of stability is¬ 
lands with the separatrices and the regions with devel¬ 
oped stochasticity. The oscillation amplitude in Figs. [4] 
(a) and (b) remains finite being of the order of the driver 
amplitude which is typical for minimal chaos [26, 27] . 

Relativistic regime 

In the relativistic regime the behaviour of the system 
0 and § is governed by two dimensionless parameters, 
a and e p . In this limit, in addition to the nonlinearity 
due to the step like profile of the restoring force on the 
l.h.s. of Eq. ([T]) , we have according to the special theory 
of relativity a nonlinear relationship between the particle 
velocity and momentum given by Eq. §. 

In the limit of small driver amplitude a <C s p the 
trajectory is confined in the vicinity of an equilibrium 
x = 0. The oscillation amplitude and frequency are of 
the order of al/2s p (in dimensional units it is eE/m e u pe ) 

and (2 Sp/l) 1 ^, respectively. In dimensional units the 
oscillation frequency is equal to uj pe . The condition of 
the displacement smallness compared to the foil width l 
is equivalent to the condition of low driver force frequency 
compared to the frequency of free oscillations. In turn it 
results in slow modulations of the oscillations. 

Fig.0 shows the case, when the driver amplitude be¬ 
ing equal to a = 0.5(7T6 p /2 ) for e p = 15.708 is close but 
lower than the resonant level. Here and below we have 
choosen the width l = 0.05, a = 15.298, and e p = 15.708. 
For uj pe /uj = 10 this corresponds to the foil width of 



FIG. 4. Driven oscillations in near-stochastic (a,b) and 
stochastic (c,d) regimes for 5 = 5 and 9, respectively, with 
l = 0.00125. a,c) The phase planes (x,p). b,d) The Poincare 
sections showing the particle positions in the phase plane 
(x,p) at the discrete time with the time step equal to the 
period of the driving force, 27r. 


the order of 10 nm. As we see, the momentum ampli¬ 
tude Pm « 0.1 is relatively small (see Figs. [ 5 ] (a) and 
(b)). It remains small for a significantly long time in¬ 
terval. The corresponding oscillation period of the order 
of 4 Pm/£p ~ 3 x 10 -2 is substantially shorter than the 
period of the driving force, which modulates the high fre¬ 
quency oscillations as seen in Fig. [ 5 ] (a). The trajectory 
fills the region (—0.003 < x < 0.003, —0.15 < p < 0.15) of 
the phase plane shown in [ 5 ] (c) due to the anharmonicity 
of the nonlinear oscillator. The driver force also mod¬ 
ulates the particle motion in the x direction as clearly 
seen in Fig. [ 5 ] (b). The Poincare section showing the 
particle positions on the phase plane (x,p) at the dis¬ 
crete time, with the time step equal to the period of the 
driving force, 27r, is plotted in Fig. [ 5 ] (d). We see that the 
trajectory is split into two sub-loops due to the oncet of 
the stochastic regime. This corresponds to a superposi¬ 
tion of the regular periodic circulation and the relatively 
weak stochastic motion. 

Under the conditions corresponding to the parameters 
e p and a above the threshold of the stochastic regime, 
a = 0.62 (tt£ p /2) the features of the driven nonlinear os¬ 
cillators are shown in Fig. [6] As one can see, at the ini¬ 
tial stage for the time interval t G [0, 60] the oscillation 
amplitude is significantly smaller than the driver ampli¬ 
tude (see also Fig. [ 7 ] below). Then it rises up abruptly. 
The maximum momentum oscillation amplitude as seen 
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FIG. 5. Driven oscillations in nonresonant case for e p = 
15.708, a = 0.52(7r£ p /2) « 13.05 and l = 0.0125. a) The 
momentum p v.s. time, b) The coordinate x v.s. time, c) 
The phase plane (x,p) for t G [0,27r]. d) The Poincare sec¬ 
tion showing the particle positions in the phase plane (x,p) 
at discrete times, with the time step equal to the period of 
the driving force, 2n. 

in Figs. [6] (a) and (b) reaches a level approximately 10 
times higher than the driver amplitude. The time de¬ 
pendence of the momentum, Fig. [6] (a), and coordinate, 
Fig. [6] (b), demonstrate the intermittency of the irregular 
behaviour. This clearly demonstrates that the intervals 
of high amplitude, low frequency oscillations with rela¬ 
tively large trajectory excursion from the x = 0 plane 
randomly alternate with low amplitude, high frequency 
motion when the trajectory is confined in the vicinity 
of x = 0 plane. The trajectory is localized within the 
domain of the phase plane shown in Fig. [6] (c). The 
trajectory tightly fills this domain demonstrating the er- 
godic property of the system under consideration. The 
domain size increases as time grows, i. e. the parti¬ 
cle momentum is gradually growing. In Fig. [6] (d) we 
plot the Poincare section showing the particle positions 
in the phase plane (x,p) at discrete times, with the time 
step equal to the period of the driving force, 27r. It has 
the form of a finite thickness web (it is also known as 
the Arnold web [28, [29]), which is broadened due to the 
stochastic nature of the particle motion. The property 
of the particle to migrate in the phase plane along the 
web is typical for the minimal chaos regimes [26]. In Fig. 
[6] (e) the Fourier spectrum of the coordinate x is pre¬ 
sented. Here the dependence of Inxn on the frequency 
£2 shows that the spectrum is continuous, which is one of 
the indications of the chaos regime m- 

The process of the particle energization in the initial 
stage is shown in Fig. [ 7 ] The intervals of regular, high 
frequency oscillations alternate with relatively short pe¬ 
riods of nonadiabatic dependence of the momentum and 



5 10 15 Q 20 

FIG. 6. Driven oscillations near the stochastic regime for s p = 
15.708 and a = 0.62(7T£ p / 2) = 15.298. a) The momentum p 
v.s. time, b) The coordinate x v.s. time, c) The phase 
plane (x,p). d) The Poincare section showing the particle 
positions in the phase plane (x,p) at discrete times, with the 
time step equal to the period of the driving force, 2tt. e) 
Fourier spectrum of the coordinate x. The dependence of 
In xq on the frequency Q is shown. 


coordinate on time. Around t = 7r, 27r, 37 t, and so on, 
the trajectory is “kicked” by the driver field resulting in 
the momentum Fig. 0(a) and coordinate]^ (b) amplitude 
increase. In the nonadiabatic phase the trajectory excur¬ 
sion along the x-coordinate is substantially larger than 
the oscillation amplitude in the adiabatic phase. Hav¬ 
ing been “kicked” several times the trajectory leaves the 
vicinity of the equilibrium plane x = 0. 

To compare the phase patterns for the parameters near 
and above the stochasticity threshold in Fig. [8] we plot 
the Poincare sections showing the particle positions in 
the phase plane (</>,p) at discrete times with the time 
step corresponding to the crossing of the x = 0 plane by 
the trajectory. In the case presented in Fig. [8] (a) for 
£ p = 15.708 and a = 0.62(77^/2) = 15.298 as in Fig. [6| 
the driven oscillator is near the threshold of the stochastic 
regime. The phase plane shows the overlapping relatively 
narrow resonances. Slightly increasing the driver ampli¬ 
tude, with £ p = 15.708 and a = 0.8(7tg p /2) = 19.74, we 
obtain the oscillations in the developed stochastic regime 











































































































FIG. 7. Initial stage of driven oscillations in the near res¬ 
onance conditions for e v — 15.708 and a — 0.62(7T£ p / 2) = 
15.298. at 0 < t < 47r. a) The momentum p v.s. time, b) 
The coordinate x v.s. time. 


as seen in in Fig. [8] (b). We see the overlapping of the 
broad stochastic layers similar to the pattern which is 
presented in Fig. 0(b)- 


Time dependence of the average kinetic energy 


In the limit of large driver amplitude, a e p , the non¬ 
linear oscillations are regular. In order to determine the 
parameters of domains where the oscillations are stochas¬ 
tic and where they are regular, we analyze the time av¬ 
erage of the kinetic energy, 


w{t m ) = 


(1 + a 2 ) 1 / 2 t 


[ tm [(i + pC0 2 ) 1/2 -i 

Jo L 


m JO 


dt. 


(48) 

Here t m is the averaging time. The electron kinetic en¬ 
ergy £ e is normalized by m e c 2 (l + a 2 ) 1 / 2 . In the case of 
particle oscillations in the driver field only, i. e. e p = 0 
which is similar to the mean electron energy introduced 
in Ref. [5], we have 


w 0 (tm ) 


7r(l + a 2 ) 1 / 2 


—E ( t m ; -a 2 ) 


(49) 


with the incomplete elliptic integral E(x\ k). In the limit 
t m —>• oc it yields for the time average of the normalized 
kinetic energy, 


w 0 


2 E(-a 2 ) 
7r(l + a 2 ) 1 / 2 


(50) 


For a> 1 we have wo ~ 2/tt. 

In Fig. [9] we plot the parametric dependences of the 
time averaged normalized kinetic energy, w. Fig. [9] (a) 
shows the function w(a) for fixed parameter e p . In Fig. [9] 
(b) we plot the function w(e p ) for fixed driver amplitude 
a. The dependences of the average normalized kinetic en¬ 
ergy on the parameters a and e p clearly demonstrate the 
stochastic regions above the threshold a = 0.62(tt£ p /2) 


with the width of the order of Sa « 10 and 


10 . 


The analysis of the time dependence of the kinetic en¬ 


ergy shows (see Fig. 10) that its growth approximately 
corresponds to the dependence oc t K . Within a relatively 



FIG. 8. The Poincare section showing the particle positions 
in the phase plane (0, p) at discrete times with the time step 
corresponding to the crossing of the x — 0 plane by the tra¬ 
jectory: a) near threshold of the stochastic regime for s v — 
15.708 and a = 0.62(7T£ p / 2) = 15.298; b) in the developed 
stochastic regime for e v — 15.708 and a = 0.8(7 T£ p /2) = 19.74. 



FIG. 9. Average normalized kinetic energy as a function of 
the driver amplitude a for the given parameter e v — 15.708, 
(a), and as a function of the parameter e p for given amplitude 
a — 15.298, (b). The averaging time is t m = 100 x 2n. 


short time interval the index k is close to 1/2 as predicted 
by Eq. (20) although there is a significant uncertainty. 
Over a relatively long time interval t E [100,1500] the av¬ 
eraged energy grows proportionally to t 1 / 3 according to 
Eq. ( |39| ) as seen in Fig. 10, where the time dependence 
of ln^/t 1 / 3 ) is plotted. 

In Fig. [9] we can see another stochastic region in the 
interval a E [50,60]. In Figs. 11 and 12 we present two 



FIG. 10. Dependence on time of the average normalized ki¬ 
netic energy logarithms for the driver amplitude a — 15.708 
parameter e v — 11.781. 
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FIG. 11. Driven oscillations near stochastic regime for e p = 
15.708 and a = 9.986 x 0.62(tt£ p /2) = 152.765. a) The phase 
plane (x,p). b) The Poincare section showing the particle 
positions in the phase plane (x,p) at discrete times, with the 
time step equal to the period of the driving force, 27r. c) 
Fourier spectrum of the coordinate x. It is shown the depen¬ 
dence of lnx(Q) on the frequency U. 


FIG. 12. Driven oscillations in regular regime for e v — 15.708 
and a — 9 x 0.62(77^/2) = 137.681. a) The phase plane (x,p). 
b) The Poincare section showing the particle positions in the 
phase plane (x,p) at discrete times, with the time step equal 
to the period of the driving force, 2tv. c) Fourier spectrum of 
the coordinate x. The dependence of lnx(Q) on the frequency 
Q is shown. 


examples of regular dynamics described by the solutions 
of Eqs. 01) for the parameters beyond the stochastic 
region, i. e. a « £ p . 


In Fig. 11 we plot the phase plane (x,p), Poincare sec¬ 
tion and the Fourier spectrum of the coordinate x for the 
driver amplitude and the parameter e p equal to 15.708 
and a = 9.986 x 0.62(7T5 p /2), respectively. The trajectory 
in the phase plane is localized in narrow width stripes. 
The Poincare section shows a stochastic layer with islands 
around the periodic trajectory. However, the Fourier 
spectrum is comprised of high order harmonics, which 
is typical for regular nonlinear oscillations. 


Slightly changing the driver amplitude to a = 9 x 
0.62(7T£ p /2) with the same parameter e p as in Fig. [Tl] re¬ 
sults in ergodic motion, when the trajectory completely 
fills the confinement domain in the phase plane (see [l2] 
(a). According to the Poincare section in Fig. 12 (b) and 
the Fourier transform, Fig. 12 (c), consisting of high or¬ 
der harmonics, the nonlinear oscillations are regular with 
no stochastic layer. 


Note that the stochastic-like regimes with developed 
stochastic layers seen in Figs. 11 (b) are not ergodic, 
i.e. the trajectory being confined within the finite size 
domain in the phase plane forms finite width stripes. It 
does not fill the domain. Contrary to this, there is an 
example when the regular trajectory fills the confinement 
domain as is seen in Fig. 12 (a). 


4. CONCLUSION 


In conclusion, we have demonstrated that a driven os¬ 
cillator with a step-like nonlinearity can exhibit regu¬ 
lar and chaotic oscillation regimes. The step-like non¬ 
linearity implies a maximal value of the restoring force. 
The threshold of the stochastic oscillation onset corre¬ 
sponds to the condition of a strong driver compared with 
the restoring force. It can be written as a relationship 
between the normalized electromagnetic field amplitude 
and the normalized surface density: a > e p . This in¬ 
equality is equivalent to the criterion of thin slab rela¬ 
tivistic transparency QME]. It also determines the pa¬ 
rameters optimal for ion acceleration by the laser light ra¬ 
diation pressure USED]- The model used above describes 
the anomalous electron heating due to the electron recir¬ 
culation around thin foil target. We note that the elec¬ 
tron recirculation during the ion acceleration by the light 
radiation pressure has been seen in theoretical work pre¬ 
sented in Refs. m and HD Eg. The electron heating 
can be enhanced for the tightly focused laser pulse |33] 
and for not perfectly aligned laser-target configurations 

m- 

Under the chaos conditions the oscillation amplitude 
grows proportionally with time to the one third power, 
which is an indication of the trajectory diffusion in phase 
space previously discussed in Ref. [8]. This result has 
been obtained by analyzing the Poincare mapping, by 
solving the KFPE diffusion equation, and by numerical 
integration of an ensemble of nonlinear oscillator equa¬ 
tions. 
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The obtained electron energy scaling given by Eq. (41), 
at the stochasticity threshold a « e p in the limit a 1 
can be approximated by 


Se 

m e c 2 


1.11 x 


"EW" 

cm 1 2 3 4 5 6 7 8 9 


1/2 

A [fJ,m](Tias [fs]) 1/3 , (51) 


where ri as is the laser pulse duration. This scaling is fa¬ 
vorable for achieving high energy ion beams generated 
under the conditions when the TNSA acceleration mech- 
anizm is realized in the laser interaction with thin foil 
targets. This shows substantially higher electron tem¬ 
perature than the scaling found in Ref. [5]. This is due 
to the electron recirculation around a thin foil target as 
has been found in PIC simulations in Refs. [35) [36] (see 
also Refs. El HI) resulting in anomalous electron heating. 
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